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We study the local classical and quantum critical properties of electron-vibration interaction, 
represented by the Yu- Anderson model. It exhibits an instability, similar to the Wentzel-Bardeen 
singularity, whose nature resembles to weakly first order quantum phase transitions at low tem- 
peratures, and crosses over to Gaussian behaviour with increasing temperature. We determine the 
dominant energy scale separating the quantum from classical criticality, study the effect of dissipa- 
tion and analyze its impact on correlation functions. Similar phenomenon should be observable in 
carbon nanotubes around local defects. 

PACS numbers: 73.43.Nq,68.35.Rh,71.38.-k,72.10.Fk 



I. INTRODUCTION 

Quantum phase transitions are intensively studied due 
to the governing fundamental physical properties and 
also because a number of highly interesting systems dis- 
play such behaviour. These occur at zero temperature 
at quantum critical points (QCP), and are dominated 
by quantum rather than classical fluctuations. Their 
study enriches our knowledge on classical and quantum 
critical behaviour, and can reveal the connections be- 
tween them in terms of quantum-to-classical mappings^. 
Quantum phase transitions were found to explain the be- 
haviour of Ge doped YbRb^Si^ and other heavy fermion 
materials^. 

The notion " local criticality" stands for quantum phase 
transitions in e.g. quantum impurity models which evolve 
only in time but are confined in space, such as the quan- 
tum phase transition in the sub-Ohmic and Ohmic spin- 
boson model 4 "^, the dissipation induced phase transition 
in a quantum bojs^ or the local quantum phase transition 
in the pseudogap Anderson models. 

Local vibrational modes due to foreign (adsorbed) 
molecules or lattice imperfections strongly influence the 
electronic transport and induce dephasing through in- 
elastic scattering 8 -. Molecular electronic device a 9 ' 10 are 
probed and controlled locally by single-molecule vibra- 
tional spectroscopy^ based on STM and inelastic elec- 
tron tunneling spectroscopy and the spectrum of molec- 
ular vibrations often indicates more complex behaviour 
than can be seen in bulk transport^, such as negative dif- 
ferential resistance and hysteresis^. Conductance mea- 
surements on mechanically controllable break junctions 
reveal the presence of local vibrational degrees of free- 
dom, when noble metal (Pt) electrodes were connected 
by a single molecule as Pt or H a 14 ! 15 . The vibrational 
mode softens after coupling it to electrons, and the vi- 
brational resonance in the conductance due these modes 
shifts to lower energies. A critical bosonic mode caused 
by the interplay of softening and dissipation contributes 
to transport down to very low temperatures. Thus, the 
detailed understanding of these modes and their local 
criticality beyond the mean field is essential^, since even 
away from the critical point, they leave their mark on the 



responses^. 

Recently the observation of strong phonon modes 
were reported in suspended carbon nanotubes (CN) by 
STMH). The electronic properties of CN are tunable by 
chemical doping or by changing the chirality of the tube, 
hence these systems are ideal candidates to study and 
control enhanced molecular vibration o 9 ' 10 , developing 
around lattice imperfections or encapsulated molecules 
via a local deformation potential. Criticality due to local 
electron-vibration interaction is regarded as the descen- 
dant of the Wentzel-Bardeen (WB) singularit y 16 ' 17 ' 18 , 
which arises in a one dimensional system of electrons, 
coupled to long wavelength phonons, which only allows 
for forward scattering. For a critical value of the electron- 
phonon coupling, the system becomes unstable and ac- 
quires a negative compressibility. The thermodynamic 
quantities and correlation functions for the electron- 
phonon system were studied near this singular poin t 19 ' 20 , 
indicating the presence of a phase transition: the diver- 
gence of the specific heat is accompanied by a collapse of 
the system induced by the strong electron-phonon inter- 
action. For a critical value of the electron-phonon cou- 
pling, the system becomes unstable and acquires a nega- 
tive compressibility. A recent study2i suggested that the 
WB singularity could be reached experimentally in thick 
carbon nanotubes due to phonons. 

Local vibrational modes and their critical properties 
are interesting for a variety of other reasons: critical 
modes are important in spintronics since they unavoid- 
ably lead to decoherence even at very low temperatures. 
Electron-vibration interactions are exceptionally impor- 
tant for molecular solids where highly energetic vibra- 
tional states greatly influence the electronic properties. 
Examples include the superconductivity with T c = 40 K 
in fullerides^ and the energetic electron-phonon side- 
bands in the excitonic excitation states of single- wall car- 
bon nanotubes 2 ^. 

This paper addresses the local criticality caused by 
electron-vibration interaction. We determine the critical 
exponents, and show that a weakly first order quantum 
phase transition governs the low temperature physics, 
similarly to certain fsing models. At high temperatures, 
it crosses over to classical, Gaussian behaviour. This 
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crossover influences the electronic properties as well e.g. 
in carbon nanotubes, which can be revealed by local spec- 
troscopical measurements-^. 



III. QUANTUM CRITIC ALITY 

Close to T 2 , the softening of the phonon frequency oc- 
curs as 



II. THE MODEL AND ITS BASIC PROPERTIES 

As schematization of local electron- vibration interac- 
tion, we start with the Yu-Anderson or single impurity 
Holstein model^i: 



'-Q 2 , (1) 



H = Y / e(k)c+c k +g d QH + {0)y{0) + ^- + 
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which describes d-dimensional electrons interacting with 
a local bosonic mode at a single impurity site with po- 
sition Q and momentum P. The model can be mapped 
onto one dimensional chiral fermions interacting with a 
single vibrational mode, and the fermionic field can be 
bosonized^^. Then, we arrive to an effective model of 
one dimensional coupled harmonic oscillators, i.e. the 
Caldeira-Leggett (CL) model^: 

OO 

H = v c f dxid^f + ^Qd.m + ^ + ^Q 2 , 
J V7r 2m 2 
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v c is the charge velocity, and g is the phase shift caused 
by ga, <fr(x) stems from the bosonic representation of the 
fermion field. This also represents the effective model for 
the large spin-boson model^I. After integrating out the 
bosonized electron field, &(x), the effective action for the 
phonon reads as 
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where w„ = 2imT is the bosonic Matsubara frequency, 
QnS are the Fourier components of Q(t). The main dif- 
ference with respect to CL is the potential renormaliza- 
tion (the — uj 2 Y /T2 term), which is avoided in CL to study 
the effect of pure dissipation. In our case, the phonon is 
expected to soften after coupling it to electrons on phys- 
ical ground, therefore such local term is present in the 
action^. As the phonon mode softens, its eigenfrequen- 
cies on the real frequency axis are given by^ 



lu p± = -ir ± ^g(i-r/r 2 )-r 2 , 



(4) 



where T 2 = ttluq/AW <C t^o <C W, W is the bandwidth of 
the conduction electrons, and T = ir(gp) 2 /2m for small 
g, and approaches r 2 as g — > oo. Here, p = 1/2itv c is the 
chiral electron density of states. The explicit dependence 
of r on ga cannot be determined by the bosonization 
approach^. The real part of the phonon frequency re- 
mains finite (underdamped) for r < T\ w r 2 (l — T^/oJq). 
For Ti < r < r 2 , the oscillatory behaviour disappears 
from the phononic response (Rea> p ± = 0), and two dis- 
tinct dampings characterize it (overdamped). 
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where y > is the distance from criticality (the effective 
reduced " temperature" ) : 
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and following Ref. Q], the characteristic energy scale is 
expected to vanish as uj p+ ~ //'. This defines v = 1, 
and using2£ v = 1/yt, the "thermal" scaling exponent is 
y t = 1. The extra i in Eq. ((5]) signals the dissipative 
nature of the transition. Note, that criticality is tuned 
by r and not by the temperature, so yt belongs to T. At 
the same time, uj p - approaches — 2ir 2 . 

Since the problem is effectively one dimensional, pos- 
sesses "0" spatial dimension and evolves only in time 
(z = I), this allows us to set d e ff = d + z = 1. The 
dynamical exponent, z can be melted in the definition 
of d, since there are no separate spatial and temporal 
dimensions. 

To proceed with the exploration of the critical proper- 
ties of our model, we evaluate its free energy. Using Eq. 
© or following Ref. HH, it is obtained as 
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where F e is the phonon free contribution of electrons. 
Focusing on the most singular contribution in y at T = 0, 
we get 
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leading to 2— a = 1 = d e ff/yt (using ln(x) = lim s ^( i (x £ — 
l)/e), from which yt = 1 is deduced in accordance with 
the exponent we got from the vanishing of the character- 
istic energy scale in Eq. ([5]) . This suggests the emergence 
of a weakly first order quantum phase transition 2 ^. Inter- 
estingly, many experimentally observed quantum phase 
transitions belong to this category^. In contrast to true 
first order transitions, where f s ~ \y\, we have logarith- 
mic corrections in y, causing the second derivative of f s 
(the equivalent of the specific heat) to diverge in a power- 
law fashion as 1/y. True first order transitions are not ac- 
companied by critical fluctuations. However, the weakly 
first order nature of the transition here allows for critical- 
ity to develop, which leaves its mark on the fluctuations. 
The mean square of the Q field follows from the effective 
action as 
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TABLE I: Summary of the critical exponent in the quantum 
and classical critical regions, 77 is only determined through 
77 = d e ff + 2 - 2y h . 



and diverges at T = as (Q 2 ) ~ |ln(y)|. The weak 
divergence of the fluctuation is a direct consequence of 
the weakly first order nature of the instability at T = 0. 
At finite temperatures, this crosses over to T/y type of 
divergence, as discussed below. 

By coupling an external field (V) to the position as 
VQ, we can study the resulting distortion within linear 
response. According to scaling 3 -!, 
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where $q(x) is a scaling function. From Eq. ((3J), we get 
V 
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from which we deduce (3 — 1/5 = 0. The scaling function 
is determined as $q(x) = 1/x. This turns out to be close 
to that of the one dimensional Ising model 3 -!. 

The divergence of the order-parameter susceptibility 
defines the 7 exponent. In the present case, the phonon 
field Q is expected to fluctuate close to the critical cou- 
pling T2, and plays the major role in the instability of 
the system. Its susceptibility is obtained as 
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with 7 = 1. This divergence is analogous to the diverging 
compressibility at the WB singularit y 19 ' 20 . It gives for 
the exponent of the external field 2 ^ 



yd + d + z 
Vh = o = !• 



(13) 



Therefore, the universality class of this problem is defined 
by d e ff = j/j = = 1, as is summarized in Table. HI 



IV. RELATION TO THE ISING MODEL 

Similar exponents characterize the one dimensional 
classical ferromagnetic Ising model as well. There, for 
any real h (longitudinal magnetic field) and finite T, 
there is no sign of criticality. However, at T = h = 0, the 
correlation length diverges, indicating the presence of a 
critical point. Following the ideas of scaling, the expo- 
nents were determined^ as /3 = l/<5 = 0, a ~ ^ = v = 
7] = 1, similarly to what we find here, and the spatial 



dimension of this classical model is d — 1. Interestingly, 
this criticality can be brought to finite temperature, if 
we consider an infinitely strong ferromagnetic chain or a 
line of defects between two neighbouring column a 32 ' 33 in 
the two dimensional Ising model on a square lattice, as 
shown in Fig. [TJ The critical behaviour on the defect 
chain differs from the bulk critical behaviour, e.g. the 
critical exponents vary with the defect chain strength. 
In the case of a chain with infinitely strong ferromag- 
netic coupling, the defect contribution to the specific heat 
diverges 3 ^ with a = 1 , and the order parameter exponent 
on the defect chain is (3 = 32 i 35 , and the other local crit- 
ical exponents also agree with that of a one dimensional 
ferromagnetic Ising chain at T — 0. Note, that the tran- 
sition temperature of the two dimensional Ising model 
is not affected by the presence of defect line. In this 
respect, the effective one-dimensional critical behaviour 
of the defect is embedded in a two dimensional critical 
region, which facilitates its observation. 



FIG. 1: (Color online) Two types of defect chains (red dashed 
and thick black lines) are shown for the 2D Ising model, lead- 
ing to the same local critical behaviour. 



CORRELATION FUNCTIONS 



The appearance of quantum criticality is further cor- 
roborated by investigating the time evolution of the cor- 
relator of the phonon field, (Q(t)Q). The mapping of our 
d-spatial dimensional quantum system (with d = 0) onto 
a d+ 1 dimensional classical one is done in the imaginary 
time path integral formalism. There, one introduces an 
extra imaginary time dimension with size 1/T, defining 
the length of the classical counterpart. Thus, as long as 
the correlation length, £ < 1/T, the system exhibits the 
previous d + 1 dimensional quantum critical behaviour. 
However, for £ > 1/T, finite size effects are important, 
quantum effect can be neglected and the system crosses 
over to d dimensional classical problem 3 ^. 

At T = 0, the decay of the correlator in imaginary time 
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is obtained from the action, Eq. ([3]) as 

oo 

cos(wr)do; 



Gq(t) = (Q(r)Q) 
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r^ 1 cos(x)dx J (yr)~ 2 for r > 2r 2 /w 2 y, 



| a; | + LU%y T /2T 2 [- ln(yr) for r < 2r 2 /w 2 ?/, 

(14) 

and Eq. (fl"4|) is valid down to r ~ l/2r2, where the u p - 
frequency starts to play its role. Scaling predicts 2 ^ that 



Gq(t) = |y| 2 (-^)/«** T {-^- t 



$ T (ry) , (15) 



where 3> r (cc) is a scaling function, and the last relation is 
obtained for our specific model. From this, we can draw 
several important conclusions. First, the scaling form 
predicts a crossover time, separating the long and short 
r regions, to scale as r* ~ 1/y ~ y~ v ', giving v = 1, in 
agreement with v = l/yt- This is in perfect agreement 
with the second integral of Eq. (|14p . where only the 
LUQyr/2T2 combination contains r and y, suggesting 
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By approaching the instability (y — > 0), r* — > oo. Sec- 
ond, scaling docs not predict additional multiplicative 
powers of y in front of the scaling function, in nice agree- 
ment with the analytical result of Eq. (TT4"]) . Third, the 
universal scaling function can be determined from Eq. 
(H3J): <f> T (x) decreases as — ln(x) for x <C 1, and decays 
algebraically as 1 jx 2 for x 3> 1 ■ The former is character- 
istic to a quantum Brownian particle^, which dominates 
the response when r* — > oo. The latter corresponds to 
the correlator of a damped harmonic oscillator. As we 
approach y — > 0, the harmonic potential flattens, and 
disappears. Then, our particle does not experience any 
confinement and performs quantum Brownian motion, 
which is not bounded, and (Q 2 ) diverges. Instead of the 
position Q, the displacement Q{t) — Q(0) keeps track of 
its dynamics^. 

The electrons also experience criticality similarly to the 
oscillator. The local charge susceptibility diverges as 1/y 
similarly to the phononic response 2 ^. The local Green's 
function of the electrons, G e (r) = (*+(0, r)^(0, 0)) de- 
cays as 1/r for long t at T = as in a local Fermi liquid, 
but for t <C t* , it changes to 
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where je ~ 0.577 is the Euler's constant. Thus, half of 
the spectral weight is lost at short times due to scatter- 
ing off the phonon close to the critical point y — 0. Con- 
sequently, the spectral function (the Fourier transform 



of G e (r)) takes half of its non- interacting value for fre- 
quencies larger than 1 /t* . This harmonizes with the be- 
haviour of the inelastic scattering rate of the electrons^, 
which reaches its maximal value at y — * for high energy 
electrons. The above results are valid for r S> 1/2^, be- 
low which the other mode (tf p _) dominates, similarly to 
Eq. (PH. 

By increasing the temperature, it is better to work in 
real times since r is restricted in a finite slab. At T > 0, 
Gq (t) decays exponentially with a correlation time given 
by £ t = l/27rT. This defines the quantum critical region, 
when besides temperature, there is no other relevant en- 
ergy scale in the problem. The diverging £ t at T —> 
agrees with the algebraic decay of Eq. (fT4"|) at T = 0. 
With increasing temperature, we cross over to the clas- 
sical critical region at 
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For T > T* . the coherence time becomes £ t = 2r 2 /u; 2 y, 
as we expect by approaching a classical transition. The 
very existence of T* proves the importance of studying 
quantum criticality. Although we might be away from the 
critical coupling, as long as the relation T < T* holds, we 
expect to observe the same quantum critical behaviour 
as in Eqs. (fT4|) and (fl~5|) . Interestingly, T* plays an 
important role at T = as well, the two regions of the 
scaling function are separated by 2ttt* = 1/T*. This 
picture is in perfect accord with critical fluctuation close 
to y = 0, which change their nature around T*. 

The knowledge of yt and yh allows us to formally deter- 
mine all critical exponents 2 ^ (see Table JJ. These expo- 
nents satisfy the scaling as well as quantum hypcrscaling 
(2 — a = (d + z)v) relations. 



VI. CLASSICAL LIMIT 

The high temperature regime of the quantum system 
can be regarded as the finite sized classical counterpart, 
with sizes smaller than the coherence length. At fi- 
nite temperatures, the singular contribution to the free- 
energy in y can be obtained from Eq. ([7]), but it is in- 
structive to follow a different approach. In the high tem- 
perature limit, the harmonic oscillator becomes classical, 
i.e. Q and P being classical variables in Eq. ([1]). Then, 
after tracing out the electronic degrees of freedom, we 
obtain the partition function for the oscillator as 
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where we integrate the exponential of the classical en- 
ergy over all of phase space (all possible momenta and 
positions), and we added a source term VQ. Then, the 
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singular part of the free-energy is 
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giving a = 2 and 7 = 1. This accounts for the instability 
in the purely classical version of Eq. ([T|) for the oscilla- 
tor, caused by the vanishing of the vibration frequency. 
Dissipation has no effect at high T, the instability is 
caused by the potential renormalization. At this point, 
we are dealing with a classical harmonic oscillator, there- 
fore d e ff = 0, and the critical theory is formally equiva- 
lent to that of a "0" dimensional Gaussian model^. The 
other exponents are determined as (3 = — 1/2 and S = — 1, 
but the v and r\ exponents are senseless, since there is no 
dimension for the spatial or imaginary time dependence. 
We mention that the criticality of the lattice version of 
Eq. |T]), the Holstein model is expected to depend on the 
dimensionality of the electrons, and belong to a different 
universality class. 



properties of a variant of the WB singularity, triggered 
by local electron- vibration interaction. The only relevant 
energy scale of the problem is identified as T* , which 
separates the quantum and classical critical regions. The 
former belongs to the universality class of weakly first 
order quantum phase transitions, while the latter is for- 
mally equivalent to a "0" dimensional Gaussian model. 
The two qualitatively different critical regimes are acces- 
sible by local vibrational spectroscopy-^, and influence 
the electronic response (e.g. dephasing^) as well. 
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Motivated by the possibility of reaching WB singular- 
ity in carbon nanotubes^i, we have studied the critical 
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